Preprint typeset in JHBP style - HYPER VERSION 



KIAS-P12036 



Non-perturbative Effect and PAMELA Limit on 
Electro-Weak Dark Matter 



Eung Jin Chun and Jong-Chul Park 

Korea Institute for Advanced Study 
Heogiro 87, Dongdaemun-gu 

Seoul 13 0-722, Korea 

Emails: |e j chunQkias . re . kr , j cparkOkias . re . kr 

Stefano Scopel 

Department of Physics, Sogang University 

1 Sinsu-dong, Mapo-gu 

Seoul 121-742, Korea 

Email: ^copeKSsogang . ac . kr 



Abstract: We discuss the non-perturbative effects on tlie anniliilation cross section of 
an Electro- Weak Dark Matter (EWDM) particle belonging to an electroweak multiplet 
when the splittings between the masses of the DM component and the other charged 
or neutral component (s) of the multiplet are treated as free parameters. Our analysis 
shows that EWDM exhibits not only the usual Sommerfeld enhancement with resonance 
peaks but also dips where the cross section is suppressed. Moreover, we have shown that 
the non-perturbative effects become important even when the EWDM mass is below the 
TeV scale, provided that some of the mass splitting are reduced to the order of a few 
MeV. This extends the possibility of observing sizeable non-perturbative effects in the dark 
matter annihilation to values of the dark matter mass significantly smaller than previously 
considered, since only electroweak-induced mass splittings larger than 100 MeV have been 
discussed in the literature so far. We have then used the available experimental data on 
the cosmic antiproton flux to constrain the EWDM parameter space. In our calculation 
of the expected signal we have included the effect of the convolution of the cross section 
with the velocity distribution of the dark matter particles in the Galaxy, showing that it 
can alter the non-perturbative effects significantly. In the case of EWDM with non-zero 
hypercharge, we have shown that the mass splitting in the Dirac dark matter fermion can 
be chosen so that the inelastic cross section of the EWDM off nuclei is allowed by present 
direct detection constraints and at the same time is within the reach of future experiments. 
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1. Introduction 

A multiplet of the Standard Model gauge group SU{2)l x U{1)y would be the simplest 
candidate of dark matter as it requires no additional 'ad hoc' interaction or couplings. Its 
electroweak interaction leads to a right thermal relic density if its mass is in the multi TeV 
range, and its stability is guaranteed automatically for a certain higher dimensional mul- 
tiplet Q . Another important feature of such an "Electro- Weak Dark Matter" (EWDM) is 
non-perturbative effects on its annihilation into gauge bosons which modifies significantly 
the tree-level results in the determination of its thermal relic density and its indirect de- 
tection rate |3| . 

In this paper, we study such a non-perturbative correction for a Majorana or Dirac 
fermion EWDM in a wide range of its mass. That is, we will consider a "non-minimal" 
EWDM allowing an unspecified non-standard cosmology for the generation of a right relic 
density and a certain discrete symmetry for its stability. In the present Universe, the dark 
matter is highly non-relativistic and thus the wave-functions of the effective non-relativistic 
two-body EWDM states can get strongly modified by the non-abelian electroweak(EW) 
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potential in the process of their pair annihilation |l|] . This is a generalization of the "Som- 
merfeld effect" well-known for a single-component dark matter carrying an abelian gauge 
charge. Recall that the non-perturbative correction enhances (reduces) the annihilation 
cross section in case of an attractive (repulsive) Coulomb potential It will be inter- 
esting to note that an EWDM exhibits not only the usual Sommerfeld enhancement and 
resonance peaks but also a suppressed cross section for a certain choice of the parameters. 
This is a realization of the "Ramsauer-Townsend (RT) effect" observed in low-energy 
electron scattering off gas atoms. We will show that such effects are caused mainly by 
the electromagnetic (EM) interaction in the two-body states of the charged components of 
an electroweak multiplet. An important feature is that the velocity distribution of dark 
matter particles has to be included in the non-perturbative calculation of the EWDM an- 
nihilation rate, as the "Sommerfeld-Ramsauer-Townsend" (SRT) resonance effect occurs 
when a certain condition is met among the model parameters, including the kinetic energy 
and so the speed of the annihilating states. 

The non-perturbative correction to the EWDM annihilation cross section is important 
to set a limit on the EWDM, the most stringent constraints coming from the PAMELA 
measurement of the cosmic anti-proton flux ||6|, ^ and from the recent FERMI-LAT mea- 
surement of the diffuse gamma-ray emission in dwarf spheroidal galaxies [^] . In the present 
paper we will focus on the PAMELA antiproton limit, deriving a constraint which is a 



slightly stronger than the FERMI-LAT limit, and comparable to the result of Ref. |10| 
with the 'MED' propagation astrophysical parameters and a fixed secondary background. 

In Section |2|, we will give a general description of various electroweak multiplet dark 
matter candidates and present formulae to calculate the non-perturbative effect, summa- 
rizing the results in Refs. ||l], ^. Depending on whether the EWDM carries a hypercharge 
or not, it can be a Dirac or a Majorana fermion. In the first case, the neutral Dirac 
components have to split into two Majorana fermions with a mass gap sufficient to sup- 
press the (inelastic) nucleonic scattering through the exchange of a Z boson. In Section 
|3|, we discuss various features of the non-perturbative correction showing the SRT effect 
for the simplest example of a triplet EWDM with no hypercharge. Section ^ discusses 
values of the mass splitting which are large enough to avoid the current direct detection 
bound but still detectable in future experiments, depending on the dark matter mass. The 
antiproton yield from the EWDM annihilation to gauge bosons is analyzed in Section ^ 
to place a limit from the current cosmic antiproton flux measurements. We then calculate 
non-perturbative annihilation cross sections for various EWDM candidates to put a mass 
limit from the PAMELA data on the antiproton flux in Section 6. We conclude in Section 





2. General EWDM and non-perturbative effect 

The dark matter particle can be the neutral component of an SU{2)l x U{1)y fermion 
multiplet. As a specific example, we will consider a vector-like (Dirac) doublet with Y = 
±1/2 (Higgsino-like) , a (Majorana) triplet with Y = (wino-like) and a vector- like (Dirac) 
triplet with Y = ±1. Note that a certain symmetry like Z2 has to be imposed for the 
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stability of these EWDM candidates. Furthermore, the dark matter component of a Dirac 
multiplet is charged under U{1)y and thus its scattering cross section with nuclei through 
Z boson exchange is far above the current limit from direct detection experiments. This 
constraint is however invalidated if there is a mass splitting in the Dirac dark matter fermion 
and thus the heavier Majorana fermion component cannot be excited by the nucleonic 
scattering of the lighter one (assumed to be the dark matter). A detailed analysis will be 
presented later to show that a mass splitting of order 0.2 MeV would be detectable while 
still allowed by the current data. Such a mass splitting can come from a higher dimensional 
operator between the EWDM and the Higgs doublet. For instance, the Higgsino-like dark 
matter multiplet, denoted by Xu = (x^^X^) a-^d Xd = iXd^Xd) chiral representation, 

allows the dimension-four operators: 

jiHuXdf, \{HdXu)\ \{HuXd){HdXu), (2.1) 

where = {H^,H^) and = eH* represent the Higgs doublets coupling to the up and 
down type quarks, respectively. Similarly, for the triplet EWDM multiplet with Y = ±1 
consisting of two chiral fermions, Xu = {xt^,xt,x'i) and Xd = (x°,Xrf,Xd~), the mass 
splitting between the Dirac pair Xu d arise from: 

^{HuHuXdf, -^{HdHdXuf, j^{HuHuXd){HdHdXu). (2.2) 

On the other hand, the wino-like EWDM multiplet, a triplet with y = denoted by 
X = (x"*" ; X*^ ) X~ ) has only one Majorana neutral component. Note that a mass splitting 
between the charged and neutral components of order 0.1 GeV arises from the electroweak 
one-loop correction Q. In the following sections, we will assume arbitrarily small mass 
gaps among the multiplet components which make a big impact on the non-perturbative 
annihilation rate of the EWDM particle. 

In the non-relativistic pair annihilation of the EWDM, the non-perturbative effect due 
to the exchange of the electroweak gauge bosons mixes together the two-body states of 
the multiplet components. In the case of the Higgsino-like EWDM, there are three states 
formed by the charged (Dirac) component and two neutral (Majorana) components: xtx'di 
XiXi and XoXo) where Xo denotes the dark matter component. For the wino-like EWDM, 
we have two two-body states: x^X~ and x^X*^- The triplet EWDM with Y = ±1 has four 
two-body states: X^^X^"' X^X^ , XiXi and XoXo- 

The Green's functions gij corresponding to the processes summarized above, where 
the indices i and j run over the two-body states of each EWDM candidate, verify the 
Schrodinger equation p|: 



1 ^^g^ 



r 



rUDM 



Qj.2 



+ Vik{r)gkj{r) = Kgij{r), (2.3) 



with moM the mass of the dark matter particle, and the boundary condition gij{0) = 5ij 
and dgij{oo)/dr = iy^mDM{K — Vii{oo))gij{oo). Here K = m/5Af/3^ is the total kinetic 
energy of the two initial dark matter particles in the annihilation process, where /? is the 
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velocity of the DM particle in the frame of the galactic halo. Then, the dark matter pair 
annihilation cross section is given by: 



avixlxl ^ AB) = 2doidl^Tf, (2.4) 

where d^j = gQj{oo) and v = 2f3 is the relative velocity between the two incident DM 
particles. Here A, B run over the gauge bosons {W~^ ,W~ , Z,^), that is, the gauge boson 
final states are AB = {W^W~ , ZZ,'jZ,^'y). Taking the normalization of the covariant 
derivative = 9^ + ig2ApT^ for each gauge boson A, the potential matrix in Eq. ( ^^ ) 
and the tree- level annihilation matrix Ty are given by |p: 

mj\r 



,2 e 

A 

where dniiQ = m^. — m^^^, and: 



ViJ{r) = 2 6m,o6iJ-a2N,NJJ2[T^j] , (2-5) 



= .n^r!.^. n^A,^B)N.N, {T^T%^ {T^T-}^„. , (2.6) 



2{1 + SAB)mjjjy^ 



, . ^^'^) / XA+XB , {xA - xb, .,, 
where j[xa,xb) = — %A/ 1 i with xa 



2 j7^2 



A 



2 16 ml^ 



Here the normalization factor Ni is 1 or y/2 for the Dirac (charged) or Majorana (neutral) 
two-body state, respectively. 



3. Sommerfeld-Ramsauer-Townsend effect 

In this section, we present a detailed study of the non-perturbative effect on the EWDM 
annihilation cross section (jy^'^+^^ = av^^ + av^^ including both final states W^W~ 
and ZZ. The wino-like EWDM system has the smallest number of states and parameters: 
two bound states {x~^X~ x^X^) mass gap between them. For this reason, 

to simplify our discussion in this section we will focus on the example of wino-like dark 
matter, unless otherwise stated. As will be discussed in Section EWDM can copiously 
produce W and Z bosons leading to a sizable contribution to the antiproton flux measured 
by cosmic-ray detection experiments such as PAMELA. This puts strong constraints on 
the masses of the EWDM particles, as will be analyzed in Section ^. 

One of the key observations in this work is that the non-perturbative effect on the 
EWDM annihilation cross section includes not only the Sommerfeld effect [Q] which in- 
duces both an overall enhancement of the cross section and resonance peaks for particular 
values of masses and couplings |]^, but also a suppression or resonance dips. We will refer to 
this as the "Ramsauer-Townsend effect" in the DM annihilation processes. The Ramsauer- 
Townsend effect is a quantum mechanical phenomenon found in the scattering of electrons 
by noble gas atoms: the collision probability reaches a minimum when the electron kinetic 
energy take a certain value Q. This is analogous to what happens to the transmission 
coefficient of a one-dimensional potential well, which is enhanced (corresponding to a van- 
ishing reflection probability) when certain conditions are met between the kinetic energy 
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of the incident particle and the potential depth and width. A similar phenomenon can 
occur in the process of EWDM annihilation in presence of a non-perturbative electroweak 
potential. In order to see how such a Ramsauer-Townsend resonance arises in this system, 
we will perform a numerical analysis of the Schrodinger equation ( p. 3D by changing the 
mass gaps, the electromagnetic, W and Z potentials, and the velocity of the DM particle. 
In particular, the latter is an important factor which can drastically change the behavior 
of the resonance peaks and dips. For this reason, in Section ^ we will include in our calcu- 
lation of the annihilation cross section a convolution over the Galactic velocity distribution 
of the impinging dark matter particles, in order to compare it to the experimental bound. 

We mention here that a plot showing dips for particular values of masses and couplings 
in the annihilation cross section of Dark Matter particles interacting in a non-perturbative 
electroweak potential was presented in Ref . . To our knowledge, this is the only instance 
where the Ramsauer-Townsend effect has been shown in the literature in the context of 
dark matter annihilation. However, the authors of Ref. |^ did not mention this effect in 
their discussion. 

3.1 Dependence on (5?ti+ 

In the determination of the non-perturbative effects of EWDM annihilation, the splitting 
between the masses of the dark matter and charged states is a crucial factor, since it controls 
the transition of the DM state to a particle able to experience the long-range effect of the 
electromagnetic(EM) interaction. The wino-like EWDM has two states and whose 
mass splitting is defined by (5m-4_ = 112^+ — m^o . In Figure ||, we present the cross sections of 
the wino-like DM annihilation to the W^W~ and ZZ final states av^^~^^^ , as a function 
of mDM for the two representative values 6m-^- = 166 and 15 MeV, where 166 MeV is the 
typical mass splitting arising from the EW one- loop correction The velocity of the 
EWDM is fixed to be the typical value of w/c = 10"'^. Our result with (5m+ = 166 MeV, 
which shows the Sommerfeld enhancement and a resonance peak within the mass range 
presented, is consistent with Ref. [|l|. For the smaller mass gap, one finds that the peak 
positions shift to smaller ttibm and a dip appears that before was missing. In this latter 
case, the smaller mass gap allows an easier access to the charged state, inducing a stronger 
non-perturbative effect which activates the Ramsauer-Townsend suppression. Actually, in 
the next subsection we will confirm that the Ramsauer-Townsend dips appear mainly due 
to the electromagnetic interaction of the charged states. 

3.2 Dependence on the EW interactions 

In this subsection, we will modify the electroweak potentials in order to see how the SRT 
effect is affected. In Figure |2|, we first plot the annihilation cross sections to W~^W~ and 
ZZ as a function of ttidm for 5m_|_ = 15 MeV with and without the EM interaction. Then in 
Figure |^, we display the annihilation cross section for the same parameters of the previous 
figure, but we vary the masses of the Z ot W boson. In particular, in Figure §(a) we take 
mz — )• n ■ mz and in Figure |^(b) mv/ — >• n ■ my/, where in both cases n = 1/3, 1, and 3. 

Figure § shows that the dip disappears when the EM interaction is turned off, con- 
firming the role played by the EM potential in the Ramsauer-Townsend effect when the 



- 5 - 




Figure 1: Annihilation cross sections of the wino-like EWDM for the mass splitting (5m+ = 
— m^o =166 MeV (blue solid) and 15 MeV (red dotted). The brown dot-dashed line shows 
the tree-level cross section. 




Figure 2: Annihilation cross sections for the wino-like EWDM with (5m+ — 15 MeV. The solid 
blue and red dotted lines show the results with and without the EM interaction, respectively. The 
brown dot-dashed line is the tree-level cross section. 

mass gap between the dark matter and charged states is small. Furthermore, one finds 
that the behavior of the annihilation cross section without the EM interaction is almost 
the same as in Figure |l|, where the influence of the EM interaction is blocked by a larger 
mass gap. This tells us that the SRT effect is insensitive to the mass gap as far as there is 
no long-range EM interaction. 

The weak interaction via the Z or W boson exchange behaves like a long range inter- 
action for a large DM mass, as can be seen in the large mass limit of Figures |l] and |2|, where 
a sizeable SRT enhancement is present. A similar effect is expected to occur if smaller Z 
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Figure 3: Annihilation cross sections for the wino-like EWDM with Sm^ — 15 MeV. Each panel 
shows the dependence on the interactions via the (a) Z and (b) W bosons. The purple dashed, 
blue solid, and red dotted lines respectively correspond to the cases of (a) n ■ mz and (b) n ■ mw 
with n = 1/3, 1, and 3. The brown dot-dashed line is the tree-level cross section. 



or W boson masses are taken. However, this may change the resonance conditions for the 
SRT peaks and dips as well. As a consequence, Figure |3| shows that the peaks move to 
smaher values of the DM mass for smaher Z and W boson masses. Furthermore, one can 
find that the SRT effect is more sensitive to a change of the W boson mass compared to 
that of the Z boson mass. This can be explained by the fact that a change of the W boson 
mass influences not only the strength of the weak interaction via the W boson propagator, 
but also the accessibility thorough an EW charged current to a state able to interact elec- 
tromagnetically. Moreover, Figure]^ (b) shows a drastic change in the Ramsauer-Townsend 
effect: in particular, the RT resonance condition cannot be met for a light W boson mass. 

From the above discussion, it is also expected that the SRT peak and dip positions 
move to lighter DM masses as the electroweak interaction strength is increased. This is 
indeed what happens in the case of the quintuplet EWDM discussed in Ref. which 
shows a RT dip at around tudm = 2 TeV for a mass gap of order 100 MeV. On the other 
hand, no RT effect for such large values of the mass gap is observed for lower-dimensional 
EWDM and a dark matter mass within the multi-TeV range. 

3.3 Amplitudes of the wave functions 

As can be seen in Figures |l], ^, and |3|, the annihilation cross sections of the EWDM show 
peaks and dips due to the SRT effect. To see their behavior in more detail, we present in 
Figure ^ the amplitudes of the wave functions doo and do-i- connecting the two bound states 
X°x° and x^X~ [see Eqs. (U) and (|2^)] of the wino-like EWDM with 5m+ = 15 MeV. 



One finds that each wave function has peaks and dips, implying both a constructive and 
a destructive resonance behavior. However, some difference is observed in the behavior of 
the peaks and dips. In particular, while peak positions coincide in the two wave function 
amplitudes, the dip positions do not. As a consequence, the dip in the annihilation cross 
section (the red-dotted line in Figure |l]) appears between the two dips in the wave functions 
(the two lines in Figure ^ and is broader. However, the dip in the cross section can be 
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Figure 4: Wave function amplitudes for the two-body states of the wino-like EWDM with (5to+ = 
15 MeV. The blue solid and red dotted lines show jdool and \do+\, respectively. 

more pronounced and even as narrow as in the case of Sommerfeld peaks in situations 
where the dips in the two wave functions are closer. 

3.4 Dependence on DM velocity 

As discussed in previous subsections, the SRT resonances occur when certain conditions 
are satisfied among the EWDM parameters such as the electroweak interaction strength, 
the mass gap and the kinetic energy (which depends on the DM velocity). So far, we 
fixed the DM velocity^ to v/c = 10~^. However, the DM particles in the halo of our 
Galaxy have a velocity distribution which is expected, for some particular values of the 
parameters, to smooth out the pattern of peaks and dips produced by the SRT effect. 
Therefore, for a real physical system it will be crucial to include the integration over 
the velocity distribution in the calculation the EWDM annihilation rate. In Figure ||, we 
present the values of terms of the relative velocity v /c of the two DM particles 

for three representative values moM = 1423, 2550, and 3380 GeV, which correspond to the 
positions of the peaks and the dips of the red dotted line in Figure [l|. As expected, the 
DM annihilation cross section shows a dependence on the velocity of the incoming DM 
particles. For the analysis of the PAMELA antiproton fiux limit in Section |^ we will use 
the annihilation cross section obtained after velocity integration according to the following 
formula: 



where Vj-ms — 

270 km/s, 

Vesc — 550 km/s and N[vesc) is the normalization constant. 



""^This is also done in other analyses, for instance, in Refs. pi, 
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Figure 5: Annihilation cross sections for the wino-like EWDM with Sm^ = 15 MeV as functions 
of the relative velocity, v/c. For each panel, the used DM mass, ttidm, corresponds to the positions 
of the peaks and the dips of the red dotted line in Figure |l|. The brown dot-dashed line is the 
tree-level cross section for each DM mass. 

3.5 Dependence on 5mjv 

Before closing this section, we will check the dependence of the SRT effect on the mass 
splitting between the neutral states. The wino-like EWDM multiplet has only one Ma- 
jorana neutral component, but, in general, EWDM multiplets can have more than one 
Majorana neutral component, whose masses are generically different from each other. In 
order to see the effect of the neutral mass splitting on the non-perturbative annihilation 
rate, we show g-y^w+zz ^.j^g Higgsino-like EWDM taking two values of the mass split- 
ting: dm^ = m^o — m^o =0.2 and 200 MeV in Figure Here the mass difference between 
the DM and charged states, 5m^ = — m^o, is fixed to be 341 MeV which is the 
typical mass splitting arising from the EW one-loop correction of Higgsino-like EWDM 
|2|. As one can see in Figure |6|, the position of the peak is shifted from moM ~ 8000 
GeV to m-DM ~ 6800 GeV when 5mN is changed from 200 MeV to 0.2 MeV. Note that in 
the previous Figure ^, the peak position moved from ttt-dm ^ 2200 GeV to mj)M ~ 1200 
GeV when (5m-(_ was changed from 166 MeV to 15 MeV. Thus, one can conclude that the 
SRT effect for EWDM is much less sensitive to the neutral mass splitting compared to the 
splitting between the DM and charged states. In the following analysis, we will fix drriN 
to 0.2 MeV (when applicable) corresponding to a situation for which direct detection is 
not excluded by the current experimental data and might yield a positive signal, as will be 
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Figure 6: Annihilation cross sections for the Higgsino-like EWDM when SniN = m^a — m^o =0.2 
MeV (blue solid) and 200 MeV (red dotted). The mass difference between the charged and DM 
states, (5m+ = — m^o , is fixed as 341 MeV. The brown dot-dashed line shows the tree-level 
cross section. 



shown in Section ^. 
4. Direct detection 

As already pointed out in Section |2|, EWDMs with non-vanishing hypercharge Y have 
large couplings to nucleons driven by the exchange of a Z boson, so their elastic cross 
sections on nuclei is excluded by Dark Matter direct-detection experiments similarly to 
what typically happens for the simplest realizations of scenarios where the Dark Matter 
particle is a massive Dirac neutrino or a sneutrino. If, however, the elastic cross section for 
the EWDM is suppressed, inelastic scattering where the DM particle makes a transition 



to a slightly heavier neutral mass eigenstate is possible |11]. Notice that the wino-like 
EWDM has a vanishing hypercharge so is not subject to constraints from direct detection. 
In the case of the EWDM with Y ^0, the neutral (Dirac) component of the multiplet is 
split into two Majorana particles (the lightest of which is the DM particle), and thus the 
elastic cross section vanishes and only an inelastic transition of the DM particle to the 
heavier mass neutral state EWDM' is allowed. In particular, among the cases discussed in 
the previous sections, this scenario can be realized for the Higgsino-like EWDM (T = 1/2, 
Y = ±1/2) or for the triplet EWDM (T = 1) with Y = ±1. 

The detection rate of inelastic scattering is suppressed by the relevant mass splitting 
6rajs[. In particular, for a given recoil energy there exists a minimal velocity for the 
dark matter Prnin below which the kinetic energy is not sufficient to allow the transition to 
the excited state: 
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In the above equation M^v is the nuclear mass and /i is the reduced mass of the DM-particle 
and nucleus. Since the incoming velocity of DM particles is bounded from above by their 
escape velocity Vesc in the Galaxy while every DM direct detection experiment is able to 
detect DM scatterings only above a given lower-energy threshold Efh, detectable values of 
the parameter dniN are bounded from above, typically a few hundreds keV depending on 
the target material and the energy threshold, while a comparison of the expected signals 
with the measured event rates allows to determine a lower bound on 5mN- Moreover, the 
EWDM-nucleon inelastic cross section by Z boson exchange is fully determined when the 
DM mass moM is fixed 

0'EWDM,nucleon^EWDM'nudeon = C — ^Y\N - {1 - 4si^)Zy , (4.2) 

zvr 

where the mass splitting between the two neutral states is neglected, c = 1 for a fermionic 
EWDM, and Z and A'^ are the number of protons and of neutrons in the target nucleus 
with mass M^v . This implies that the allowed range for dm jy can be plotted as a function of 
rriDM- This is done in Figure |^ for the case of the latest constraints from the XENONIOO 



experiment |12]. In this figure, the solid lower curve refers to the case Y = ±1 while the 
dashed lower curve corresponds to y = ±1/2. As explained above, values of Smjq above 
the upper curve are non-detectable, since in this case cPmin > Vmax = Vesc + i^earth-"^ This 
is a purely kinematic constraint that does not depend on the EWDM-nucleus cross section, 
so the upper solid curve in the figure is common to the cases with Y = ±1 and Y = ±1/2. 
On the other hand, the lower curves represent the 90% C.L. lower bounds on (5mjv obtained 
by applying Yellin's maximal gap procedure [^] to the spectrum (consisting of two nuclear 
recoil candidates at recoil energies 7.1 keV and 7.8 keV) detected by XENONIOO in the 



range 6.6 keV< E^- < 43.3 [12|. In the calculation, we have assumed the standard value 
/Odm=0.3 GeV/cm'^ for the DM density in the neighborhood of the Sun and a Maxwellian 
velocity distribution truncated at Vesc=550 km/s with Vrms='2'^0 km/s.^ 

As it will be shown in detail, antiproton fluxes are almost insensitive to the particular 
choice of SniN, while are very sensitive to the other mass splittings. Nevertheless, in the 
following sections we will fix 6mi\f =200 keV (when applicable) corresponding to a situation 
for which direct detection is not excluded by present constraints and might yield a positive 
signal. 



5. Constraints from antiproton fluxes 

Since the EWDM is SU{2) charged, their annihilations are expected to produce W/Z 
bosons copiously whenever this channel is kinematically allowed, leading to a sizeable 

■^Notice that Pmin is defined in the detector rest frame, whiie Vesc is defined in the Galactic rest frame. 
For the Galilean boost we have assumed Vsarth~'232 km/s. 

^'In Ref. the DM Region Of Interest (ROE) 6.6 keV< Er < 43.3 is used when analyzing the data 
with a Profile Likelihood method, but is reduced to 6.6 keV< Er < 30.5 keV when applying the maximum- 
gap method. While this does not imply a significant change in the limit for elastic scattering, the inelastic 
scattering bound is very sensitive to the upper bound of the ROI. For this reason, we derive our limit using 
the whole energy range of the XENONIOO measurement. 
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Figure 7: Values for the mass splitting Sm^ which are detectable and experimentally allowed by 
the XENONIOO direct detection experiment as a function of the DM mass rriBM- The solid 
lower curve refers to the case Y = \ and the dashed lower curve to the case Y = 1/2, while the solid 
upper does not depend on the value of the cross section and is common to both cases. When 5miq 
is above the upper curve the recoil energy is always below the XENONIOO lowest-energy threshold 
Eth=^ keV. On the other hand, values below the lower curves are excluded at 90% C.L. by the 2 
nuclear-recoil candidate events observed by XENONIOO in the range 6.6 keV< < 43.3 keV by 
applying Yellin's maximal gap method p^ . 



primary contribution to the antiproton flux detected by experiments measuring cosmic- 
rays. The antiproton primary contribution from DM annihilation must be summed to 
the secondary antiproton contribution produced by energetic cosmic rays impinging on 
the interstellar medium. Although still affected by uncertainties, the latter contribution 
can be calculated in a relatively reliable way, and is in fair agreement to observation. This 
implies that no much room is left for the additional contribution from DM annihilation, and 
antiproton data can put constraints on the EWDM, namely on the EWDM annihilation 
cross section-times-velocity av. 

This is shown in Figure where the circles represent the top-of-atmosphere antiproton 
flux as measured by PAMELA , while the dashed line is the secondary flux as calculated 
in Ref. and rescaled by an overall factor 0.84. In this way, the model fits the data 
particularly well, with = 12.1 with 23 degrees of freedom. For this reason, we adopt 
this model as an estimation of the secondary flux. As far as the primary antiproton flux 
is concerned, we have used both the antiproton yields per annihilation and the propaga- 
tion model according to Ref. |15|, adopting for the latter an Einasto density profile with 
median values of the propagation parameters. Notice that, since the antiproton yields 
per annihilation corresponding to the two different final states WW and ZZ are practi- 
cally undistinguishable, in the calculation of the expected signal the total cross section 
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Figure 8: The circles show the antiproton top-of-atmosphere flux measured by PAMELA |j as 
a function of the antiproton kinetic energy. The dashed line represents the secondary flux as 



calculated in |14 and rescaled by an overall factor 0.84. The solid lines show three expected fluxes 
from DM annihilation calculated as described in Section ^ for TO£ijv/=200 GeV, 500 GeV, and 1 TeV, 
all corresponding to x^=44.2. 

^WW _|_ ^zz -g factorized. 

The result of our analysis is shown in Figures |^ and ^. In particular, for each choice of 
rriDM and av^'^ + av^^ we sum the primary and secondary contributions of the expected 
antiproton flux and use the PAMELA data points to calculate a x^: fo^^ which we assume 
an upper bound <44.2, corresponding to the 99.5% C.L. with 23 degrees of freedom. In 
this way, we obtain the solid line shown in Figure ^. In Figure ^, the three solid lines show 
the expected antiproton flux for m£)M=200 GeV, 500 GeV and 1 TeV when av^^ + av^^ 
lies on the boundary given in Figure ^. 

6. PAMELA limit on various EWDM 

As shown in Section ^, W and Z bosons produced by the annihilations of EWDMs provide 
a sizable contribution to the cosmic ray antiproton flux measured by PAMELA. Thus, 
for various EWDM models we will examine the DM mass ranges satisfying the PAMELA 
antiproton flux bound on as obtained in Section ^. As explained in Section 

3.4, the annihilation cross section of the EWDM shows a DM velocity dependence. In this 
section, we will therefore calculate the annihilation cross section considering the velocity 
integration effect, and also, for comparison, present the result with a fixed velocity, v/c = 
10-^. 



6.1 Higgsino-like EWDM 

We first show the results for the Higgsino-like (doublet with Y = ±1/2) EWDM which is 
the smallest multiplet. In Figure IC, we present the annihilation cross sections (blue solid 
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Figure 9: Estimation of the 99.5 % C.L. upper bound on the annihilation cross section times 
velocity (av) as a function of the DM mass niDM- The values of (av) above the solid line have 
>44.2 with 23 degrees of freedom when compared to the PAMELA data 0. 



lines) for two representative values 5m-^ = 341 and 8 MeV as a functions of jtidm, where 
the velocity integration effect is included. For comparison, we also show our results for a 
fixed relative DM velocity v/c = (red dotted lines). As already explained, in this plot 
and in the following ones the neutral mass splitting Smi\i is taken to be 0.2 MeV, when 
applicable. For the higher mass splitting Sm^ = 341 MeV (a), we extended the DM mass 
range up to 10 TeV to see the first Sommerfeld peak. If we further extended the mass 
range, Ramsauer-Townsend dips would appear as well. On the other hand, for the smaller 
mass gap (5m_|- = 8 MeV (b), the SRT resonances appear in the smaller DM mass range as 
discussed in Section 3. 

As shown in the figures, the velocity integration smooths out the peaks and dips and 
changes the positions of the peaks in (b). With only the tree-level cross section, the region 
"T-DM ^ 364 GeV is ruled out by the current PAMELA data. However, the annihilation 
cross section is enhanced by the SRT and consequently the excluded region is a bit extended 
to mDM ^ 382 GeV in the case of the typical charged mass splitting for the Higgsino-like 
EWDM, (5m+ = 341 MeV ||. In the case of a smaller mass splitting, the PAMELA Hmit 
can constrain also small bands around the Sommerfeld peaks, in addition to the low mass 
region, as can be seen in Figure |l^ (b). 

6.2 Wino-like EWDM 



Figure 11 shows the annihilation cross sections for the wino-like (triplet with Y = 0) 
EWDM in the same way as in the Higgsino-like EWDM case. The two values 6m^ = 166 
and 6 MeV are taken for the analysis, in order to show the dependence on the charged mass 



splitting. One can see from Figure 11 that the velocity integration makes a big change in 




Figure 10: Annihilation cross sections to W^W~ and ZZ for the Higgsino-like EWDM with 
dm+ — (a) 341 MeV and (b) 8 MeV, and SniN — 0.2 MeV. In each panel, the blue solid line is our 
final result obtained after velocity integration, while the red dotted line is the cross section with 
the fixed velocity, v/c= 10~^. The black long-dashed line shows the upper limit obtained from the 
PAMELA antiproton fiux data analysis in Section ||. The brown dot-dashed line is the tree-level 
cross section. 
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Figure 11: Annihilation cross sections to W^W and ZZ for the wino-like EWDM when (5m+ = 
(a) 166 MeV and (b) 6 MeV. Each line represents the same thing as in Figure |l^. 

the case of the smaller mass gap, erasing out some of the peaks and dips in ^ (b). While 
the tree-level result excludes the mass range w-dm ^ 533 GeV, the non-perturbative effect 
extends the constrained region up to w-dm ~ 664 GeV for the representative charged mass 
splitting of the wino-like EWDM, (5m+ = 166 MeV[§]. For the smaller mass gap 5m+ = 6 
MeV the PAMELA limit gets stronger, excluding DM masses below 900 GeV and also 
bands around the peaks which are larger compared to the Higgsino-like EWDM case. 

6.3 Hyper-charged Triplet EWDM 

Now let us consider a more complicated case, the triplet EWDM with y = ±1 which 
has one doubly charged, one singly charged, and two neutral components. Thus, the DM 
component Xo can have three mass splittings: (5m++ = m^++— m^o, (5m_(_ = m^o, and 

(JmAT = m^o — JTi^o. Figure 12 presents the annihilation cross sections for (5m++,(5m+) = 
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Figure 12: Annihilation cross sections to W^W~' and ZZ for the triplet EWDM with Y = ±1. 
The used values of ^m^+ ^^Ar are shown in each panel. Each line represents the same thing as in 
Figure |o[ 

(1400, 525), (100, 525), (1400, 15), and (100,15) MeV fixing 6mN = 0.2 MeV. Note that 
(5m++ = 1400 MeV and 5m+ = 525 MeV correspond to the typical mass splittings due to 
the EW one- loop corrections [||. 

One can see the dependence of the SRT effect on the doubly charged mass splitting 
5m_|.+ comparing the two cases (5m++, (5m4_) = (1400,525) and (100,525) MeV in Figures 
[l^ (a) and (b), while the dependence on the singly charged mass splitting 5m+ can be 
understood comparing the two cases ((5m++,5m+) = (1400,525) and (1400,15) MeV in 
Figures |l^ (a) and (c). In particular, one can see that a larger number of peaks and dips 
appears when 6m^-^- is reduced by a factor 14 compared to the case when 6m^ is reduced 
by a factor 35. This shows that the SRT effect is more sensitive to (5?n++ than to 5m+ 
due to the stronger EM interaction of multiply-charged states. In addition. Figure ^ (d) 
shows the combined effect of the changes of 5m++ and 6m^ shown separately in Figures 
[l^ (b) and (c). The hyper-charged triplet EWDM has stronger EW interactions compared 
to the Higgsino-like or wino-like EWDM and thus exhibits a stronger SRT effect. As a 
consequence, the excluded mass region reaches about 3 TeV for the typical mass gaps of 
(5m++ = 1400 MeV and 5m+ = 525 MeV, while the tree-level limit is as low as about 800 
GeV as shown in Figure |l^ (a) . It is interesting to see that regions of small DM mass are 
allowed for smaller mass gaps, as shown in Figures |l2| (b,c,d) due to the RT effect. One 
would expect to find more regions of lower DM mass allowed for lower mass gaps. 
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7. Conclusions 



In the present paper, we have discussed the non-perturbative effects occurring in the anni- 
hilation cross section of an "Electro- Weak Dark Matter" (EWDM) particle belonging to an 
SU{2)l X U{1)y multiplet, when the splittings between the DM state mass and that of the 
other charged or neutral component (s) of the multiplet are treated as free parameters. In 
particular, we have considered a vector- like (Dirac) doublet with Y = ±1/2 (Higgsino-like) , 
a (Majorana) triplet with Y = (wino-like) and a vector- like (Dirac) triplet with Y = ±1. 
In all these examples, an ad hoc symmetry has to be imposed for the stability of EWDM, 
and we have allowed for an unspecified non-standard cosmology for the generation of the 
right dark matter relic density, since the thermal abundance of EWDM is typically below 
that required by observation unless its mass is in the multi-TeV range. Moreover, in the 
case of EWDM charged under U{1)y, severe constraints from direct detection searches 
exist on the elastic cross section off nuclei. However, these limits can be circumvented in 
presence of a sufficiently large mass splitting (of the order of 0.2 MeV) in the Dirac dark 
matter fermion, so that only inelastic scattering is allowed and kinematically suppressed. 

As a result of our analysis, it is shown that EWDM exhibits not only the usual Som- 
merfcld enhancement of the cross section with resonance peaks at particular values of the 
dark matter mass, but also a suppressed cross section for particular choices of the parame- 
ters. The latter phenomenon is a realization of the "Ramsauer-Townsend effect" observed 
in low-energy electron scattering off gas atoms. Moreover, we have shown that the EWDM 
mass for which non-perturbative effects become important is particularly sensitive to the 
mass splittings between the dark matter and the charged components of the EW multiplet, 
and is driven below the TeV scale when this mass splitting is reduced to a few MeV. In 
particular, we have shown that when the mass splitting gets smaller the transition of the 
dark matter particles to electrically charged states is made easier, and it is the electromag- 
netic long-range interaction between these charged states which is responsible both of the 
Sommerfeld enhancement of the cross section and of the Ramsauer-Townsend suppression, 
even when the dark matter mass is not much larger than the EW gauge boson masses. No- 
tice that only mass splittings larger than 100 MeV, induced by EW radiative corrections, 
have been considered so far in the literature, so that, before our analysis, the Sommerfeld 
effect had been discussed only in the context of multi-TeV scale dark matter. 

Based on the results explained above, we have then used available experimental con- 
straints on the exotic component of the antiproton flux in cosmic rays to put constraints 
to the EWDM parameter space. Since non-perturbative effects depend on the velocity of 
the dark matter particles, we have calculated the annihilation cross section considering the 
effect of the convolution of the cross section with the velocity distribution of the dark mat- 
ter particles in the Galaxy, showing that in some cases this can smear out or significantly 
modify the pattern of Sommerfeld-Ramsauer-Townsend peaks and dips obtained for a fixed 
value of the velocity. Typically, wc have found constraints on the EWDM mass ranging 
between a few hundred GeV to a few TeV, depending on the specific EWDM realization 
and on the values of the mass splittings. However, wc have also found that, by an ap- 
propriate choice of the mass splittings, sometimes the Ramsauer-Townsend suppression in 



-17- 



the annihilation cross section can ahow to recover narrow intervals at lower values of the 
EWDM mass. 

Finally, in the case of EWDM charged under U{1)y we have found that the phe- 
nomenology described above is not particularly sensitive to the mass splitting between the 
two neutral Majorana states. As a consequence, this mass difference can be chosen so that 
the inelastic cross section of the EWDM off nuclei is allowed by present direct detection 
constraints and at the same time is within the reach of future experiments. 
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